Parity Violation in Proton-Proton Scattering at 221 MeV 
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TRIUMF experiment 497 has measured the parity violating longitudinal analyzing power, A z , 
in pp elastic scattering at 221.3 MeV incident proton energy. This comprehensive paper includes 
details of the corrections, some of magnitude comparable to A z itself, required to arrive at the final 
result. The largest correction was for the effects of first moments of transverse polarization. The 
addition of the result, A z = (0.84 ± 0.29(stat.) ± 0A7(syst.)) x 10~ 7 , to the pp parity violation 
experimental data base greatly improves the experimental constraints on the weak meson-nucleon 
coupling constants h p p p and hPJ" , and also has implications for the interpretation of electron parity 
violation experiments. 

PACS numbers: 21.30.Cb, ll.30.Er, 24.70.+S 25.40.Cm 



I. INTRODUCTION 



This experiment determines the parity-violating longi- 
tudinal analyzing power, A z = (a + — o~)/(o + +<7~), in 
pp elastic scattering, where er + and a~ are the scatter- 
ing cross sections for positive and negative hclicity. The 
measurements were performed in transmission geometry, 
with beam energy and detector geometries selected to en- 
sure that parity mixing in the lowest order 1 Sq — 3 Pq par- 
tial wave amplitude did not contribute to the measured 
A z , hence leaving a parity violating asymmetry arising 
almost entirely from 3 P-2 — 1 D2 mixing [[J. This ampli- 
tude has not previously been studied experimentally, and 
the possibility is unique to the energy regime accessible 
with the TRIUMF cyclotron. The energy at which the 
contribution to A z from the lowest order 1 S'o— 3 Po mixing 
vanishes is determined by the well known strong nuclear 
phase shifts. The scale factors multiplying this and other 
partial wave contributions are set by the weak interac- 
tion, pp parity violation experiments determine these 
scale factors experimentally. In the context of the weak 
meson exchange model Q , the TRIUMF measurement of 
A z determines primarily the weak p-meson-nucleon cou- 
pling constant h pp = (h® + /i 1 + hZ/VE), where the su- 
perscripts refer to isospin change [3J. Precision results 
already obtained by the SIN group at 45 MeV 3 and 
the Bonn group at 13.6 MeV [5| determined essentially 
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FIG. 1: Principle of the measurement. A longitudinally po- 
larized proton beam was passed through a liquid hydrogen 
target. The beam current before and after target was de- 
tected and the signals were subtracted. The presence of a 
component of the difference signal synchronized with spin flip 
indicates a parity violating dependence of the transmission 
through the target on the helicity of the incident protons. 



the sum hf + h p J, where = h% + h l u . With the ad- 
dition of the TRIUMF result at 221.3 MeV, h™> and h p f 
are now determined separately for the first time. 

It is important to have experimentally determined val- 
ues of the weak meson-nucleon couplings, as theoreti- 
cal calculations of their values are quite uncertain and 
the correct values are helpful in the interpretation of 
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FIG. 2: General layout of the TRIUMF parity experiment. 
(OPPIS: Optically Pumped Polarized Ion Source; SOL: Spin 
Precession Solenoid; FCSM: Ferrite Cored Steering Magnet; 
IPM: Intensity Profile Monitor; PPM: Polarization Profile 
Monitor; TRIC: Transverse Field Ionization Chamber) 

the results of other parity violation experiments. They 
are needed, for example, in calculations of the proton 
anapole moment [(|, one of the radiative corrections to 
the electron-nucleon isovector axial form factor in exper- 
iments such as SAMPLE and G° 0. 



II. BACKGROUND 

TRIUMF E497 @, was first funded in 1988, and the 
required new beamline was completed in 1994. A major 
effort to understand and minimize systematic error con- 
tributions was then undertaken. Following many years of 
effort that resulted in the reduction of both the helicity 
correlated beam modulations Axi and the sensitivities 
to these modulations, the first significant data set 
for E497 was acquired in February and March of 1997, 
with a statistical error on A z of ±0.5 x 1CP 7 , and most 
systematic errors at or below the 10~ 7 level. That result 
represented a major milestone for the experiment jlCj . 
Data taking continued in 1998 and 1999, the final re- 
analysis of all the data was completed in early 2001 [Tlj . 
and the result for A z was published ^3- The present 
paper presents detailed descriptions of the experiment, 
the data analysis, and systematic error corrections, that 
were not included in the Letter IFJ. 



III. EXPERIMENTAL SETUP 
A. General 

The principle of the experiment is straightforward. A 
longitudinally polarized proton beam is passed through 
a liquid hydrogen target and the change in transmission 
when the spin of the incident protons is reversed is the 
parity violating signal. To measure this, the beam cur- 



Vert. 1 gate 
Hor. 1 gate 
Vert. 2 gate 
Hor. 2 gate 
TRIC/IPM gate 





1 


1 

2 1 3 
1 


4 


5 


1 

6 1 7 
1 


8 




+ 


1 


+ 


+ + 


























I I U I 

I I I 1 


1 II 1 1 




p p 1 p 1 1 p 1 




n a 1 1 1 1 1 1 1- 




-r- b b 1 ~^^^^r- 




— i — : — i — iu p : — 1 — r 




ill I U U i 







TIME (ms) 25 50 75 100 125 150 175 200 

FIG. 3: Sequence of events in one 200 ms, 8-state cycle. Each 
state starts with a PPM scan. In the first two states of an 
octet, PPM1 scans vertically, on the next two states, PPM1 
scans horizontally. This 4-state sequence is then repeated 
with PPM2. 



rent before and after the liquid hydrogen target was mea- 
sured and the signals were subtracted (Fig. HJ. The 
spin state was then changed in a special pattern at 40 
spin states per second and synchronous detection was 
used to extract that component of the difference signal 
that was synchronized with the spin flip. Unfortunately, 
beam parameters other than spin changed when the spin 
was flipped, and great pains had to be taken to measure, 
and correct for, the false signals resulting from these un- 
wanted helicity correlated changes. Technical details of 
the systems required to do this will be described in de- 
tail in what follows, but an idea of the complexity can 
be obtained from Fig. [21 identifying the major pieces of 
equipment, and Fig. showing the sequencing of the 
control measurements within a 200 ms, eight spin-state 
cycle. 

Because of the importance of understanding systematic 
errors in such an experiment, essentially all the pieces of 
equipment, from the ion source to the beam dump, had 
to be considered integral parts of the parity experiment. 
Figure [2 shows the main subsystems of the experiment - 
the TRIUMF optically pumped polarized ion source (OP- 
PIS), the cyclotron, the beam transport, and the special- 
ized parity instrumentation . A 5 ^A transversely po- 
larized beam was transported to the cyclotron through an 
approximately 50 m long injection beamline. A 200 nA 
beam at 75 - 80% vertical polarization was extracted at 
221.3 MeV. Spin precession through pairs of solenoid and 
dipole magnets resulted in delivery of a longitudinally 
polarized beam to the parity apparatus. There were two 
complementary states of the spin transport, the "positive 
helicity" and "negative helicity" beamline tunes, which 
transported spin-up in the cyclotron into either + or — 
helicity at the parity apparatus. 

In the last section of the beam line, the longitudinally 
polarized beam first passed through a series of diagnostic 
devices - two beam intensity profile monitors (IPMs) and 
two transverse polarization profile monitors (PPMs) - 
before reaching the LH2 target, which was preceded and 
followed by transverse electric field, parallel plate ioniza- 
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tion chambers (TRICs) to measure the target transmis- 
sion. A third IPM was located immediately in front of 
the LH2 target, inside the cryostat vacuum. 



B. Optically Pumped Ion Source 

In the ion source, 9 W of 795 nm laser light was used to 
optically pump a rubidium vapor whose polarization was 
ultimately transferred to the protons of the H~ beam. 
The polarization was reversed using small tilting ctalons 
to make rapid frequency adjustments to the two pump- 
ing lasers to match either the a + or o~ component of 
the rubidium Di transition; the two are separated by 
93.5 GHz. No macroscopic electric or magnetic fields 
were altered. This minimized helicity correlated changes 
in accelerated beam parameters other than polarization. 
This was very important to the success of the experi- 
ment. The ion source is the ultimate origin of all helic- 
ity correlated modulations, and the more that unwanted 
modulations could be reduced at the source the less cor- 
rections were required later. The extracted beam current 
required by the experiment was not large (200 nA) so, to 
reduce unwanted helicity correlated modulations, it was 
possible to sacrifice most of the OPPIS intensity in re- 
turn for beam quality. For example, the RF bunchers 
in the injection beamline, which enhance the cyclotron 
transmission by a factor of five, also amplify the coher- 
ent energy modulations by two orders of magnitude, and 
could not be used during the parit y exp eriment. High- 
current OPPIS development |15l iWL Il7j for high energy 
accelerators proceeded at the same time as development 
of the TRIUMF source and contributed greatly to the 
parity violation experiment. 

The polarization of the rubidium vapor was monitored 
and controlled on-line by observing the Faraday rota- 
tion of light from an additional 100 mW, TiS probe laser 
that emitted 780.8 nm light - close to the D2 transition 
of rubidium. The polarization of the linearly polarized 
probe laser light rotated through an angle proportional 
to the rubidium vapor polarization. The Faraday mea- 
surements also provided confirmation of the helicity state 
of OPPIS. The Faraday rotation signal was encoded as a 
frequency to prevent helicity correlated signals from be- 
ing present in the electronics racks. Details of the OPPIS 
Faraday rotation system are described elsewhere 0] ■ 

Every effort was made to tunc OPPIS for minimum 
intensity change on spin flip. The main technique used 
to do this was to keep the rubidium polarization, as mea- 
sured by the Faraday rotation, at close to 100% and 
with the two spin states matched to better than ±0.5%. 
It was not possible to eliminate helicity correlated cur- 
rent modulation completely, but under normal data tak- 
ing conditions, helicity-correlated current modulations 
AI/I = (1+ - + I~) of a few parts in 10 5 were 

routinely achieved. 

To enable measurement of the sensitivity of the exper- 
iment to coherent intensity modulation, provision was 



made at the ion source to produce an intentional, con- 
trolled, intensity modulation when desired. This was 
done using an auxiliary 18 W green argon laser (Spectra 
Physics Inc., operating on all lines with most power at 
514 nm and 488 nm) beam that co-propagated with the 
H ~ beam along the 30 m long horizontal section of the in- 
jection beam line, neutralizing through photodetachment 
a small fraction of the beam. The photodetachment laser 
could be interrupted synchronously with the parity spin 
sequence, so that the beam current in every other "spin 
off" (i.e. with the optical pumping lasers blocked with 
a shutter) data taking cycle was modulated at the 0.1% 
level. 



C. Cyclotron and Beamline 

1. Beam Transport 

After extraction from OPPIS, the H~ beam passed 
through a Wicn filter that was tuned to produce vertical 
polarization at the entrance to the cyclotron. The injec- 
tion line from the ion source to the cyclotron used electro- 
static elements and was magnetically shielded from the 
fringe field of the cyclotron. Instabilities in the polar- 
ized ion source, beamline power supplies, and mechan- 
ical vibrations of the whole injector building structure, 
cause the beam position at the injection point to fluctu- 
ate. These fluctuations arc converted to energy and cur- 
rent modulations of the accelerated beam. These were 
significantly reduced by a position stabilization feedback 
system installed in the vertical section of the injection 
beamline 0. The system was based on two split-plate 
beam position monitors, with correction voltages applied 
to electrostatic steering plates. About 50% of the beam 
was lost on the split-plates, but the result was a signifi- 
cant improvement in the beam stability, and a reduction 
in noise. The sampling rate of the integrated current 
feedback amplifier was 1 kHz, so spin-flip correlated po- 
sition modulations produced in the source at 40 Hz were 
also reduced by the position stabilization system. 

After injection into the cyclotron, the beam was ac- 
celerated to 221.3 McV in the cyclotron, and was ex- 
tracted by a thin stripping foil. Various stripping foil 
designs were tried, but most of the data were taken with 
a 2.5 mm wide x 26 mm high x 5 mg/cm 2 thick, py- 
rolytic graphite foil. This foil was mounted in a special 
"bow-saw" shaped holder that supported the foil from 
both ends to prevent curling. Following extraction from 
the cyclotron, the proton polarization vector was pro- 
cessed through ± ~ 63° in the first solenoid, 88.61° in 
the first (40°) dipole, ± ~ 87° in the second solenoid, and 
26.58° in the final (12°) dipole, resulting in a longitudi- 
nally polarized beam that was transported to the parity 
apparatus. The sign of the solenoid rotation was cho- 
sen depending on the desired helicity of the tune and the 
exact solenoid strengths were fine tuned empirically to 
produce pure longitudinal polarization in the presence of 
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TABLE I: Beam energy measurements made during the three 
running periods. The uncertainties shown are statistical only. 
The absolute calibration of the BEM is known to ±0.23 MeV. 



Running period 


beam energy 


1997 


221.30 ± 0.03 MeV 


1998 


221.5 ±0.1 MeV 


1999 


221.3 ±0.1 MeV 



the cyclotron fringe field. In contrast with systems using 
only one solenoid and one dipole, which can produce lon- 
gitudinal polarization only at one energy, the TRIUMF 
system is capable of producing longitudinal polarization 
at any proton energy up to 500 MeV. [l9j j. 

Ideally, the beam transport should produce an achro- 
matic waist downstream of the LH2 target, to minimize 
the effect of first moments of transverse polarization, 
which are the dominant source of systematic error in this 
experiment, as they couple to the large (~ 0.3) parity 
allowed transverse analyzing power A y . The rotation of 
phase space introduced by the superconducting solenoid 
magnets was approximately compensated for by rotating 
the quadrupole doublet before the final 12 degree bend 
by 37.4° and the quadrupole triplet following the bend 
by 17.2°. To reverse the spin direction at the parity tar- 
get relative to the spin in the cyclotron, the directions 
of the solenoidal fields and the signs of the quadrupole 
rotations were reversed. In principle, that should have 
been be all that was required. In practice, small adjust- 
ments to some quadrupole and steering magnets were 
necessary, but empirical tunes were developed that kept 
these changes very small (1-2%). 

Between the second solenoid magnet and the final 12° 
bend, where the polarization of the beam has both a lon- 
gitudinal and a sideways component, a four branch po- 
larimeter [2(j was used to measure the transverse polar- 
ization components at regular intervals. Since the angle 
of the polarization vector is known at this location, the 
absolute beam polarization could be determined. Once 
the beam passed the final dipole magnet, the polariza- 
tion was longitudinal and the PPMs determined the small 
transverse components. The upstream polarimcter was 
also used, with purely vertically polarized beam, to check 
the absolute calibration of the PPMs. 

Just after the upstream polarimeter, a Beam Energy 
Monitor (BEM) [2l| was available to measure the en- 
ergy of the proton beam. The BEM achieves a statistical 
precision of ± 20 keV in approximately one hour. The 
absolute calibration is known to ±230 keV j^. The 
(BEM) target could not be inserted during normal data 
taking, so periodic BEM measures were made to check 
the beam energy. Table [I] summarizes the beam energy 
measurements. 



2. Beam Energy 

The beam energy of 221.3 MeV at the target entrance 
was chosen so that the ( 1 S'o — 3 Pq) contribution would 
be zero when the acceptance of the detectors and the 
energy loss in the target was taken into account. This 
was done using a Monte Carlo program that simulated 
the experiment using detailed target and detector geome- 
tries, and including the energy dependence of the strong 
pp phase shifts. It was found that the beam energy was 
not overly critical; a departure of one MeV from the en- 
ergy for which the ^Sq — 3 Pq) contribution integrated 
to zero caused this mixing to contribute only ~ 10 -9 to 
A z . Theoretical calculations of A z based on total cross 
section, and which assume no energy loss in the target, 
show the ( x 5 - 3 P ) zero at 225 MeV, where the ^0 
and 3 Po strong phase shifts are equal and opposite [23|. 
According to the Monte Carlo simulation, this theoreti- 
cal A z at 225 MeV should be compared to the TRIUMF 
A z multiplied by 1.02 ± 0.02. 

D. Specialized Instrumentation 

1 . General 

Very strict constraints had to be imposed on the inci- 
dent longitudinally polarized beam to limit both random 
and helicity correlated modulations of the beam inten- 
sity, energy, horizontal (2) and vertical (y) position and 
direction, beam width (cr x and a y ), transverse polariza- 
tion (P x and P y ), and first moments of transverse polar- 
ization ((xP y ) and (yP x )). The approach followed was 
to design the experimental apparatus for minimum sen- 
sitivity to beam property variations, to monitor helicity- 
correlated beam properties during data taking, and to 
correct the data for all significant helicity correlated ef- 
fects. A program of auxiliary calibration measurements 
was interwoven in the regular data taking cycle in order 
to establish the sensitivity of the apparatus to all mea- 
surable systematic effects under data taking conditions. 

The custom-built parity instrumentation occupied ap- 
proximately the last 10 m of the bcamlinc, between the 
last quadrupole magnet and the west wall in the TRI- 
UMF proton hall (beamline 4A/2). Transverse field par- 
allel plate ion chambers TRIC1 and TRIC2 measured 
the beam current incident on and transmitted through 
the target. The parity violation signal was derived from 
the helicity-corrclatcd analog signal difference between 
the beam currents measured by the two TRICs. Up- 
stream of the target were two polarization profile mon- 
itors (PPMs) to measure the distributions of transverse 
polarization P y (x) and P x {y) across the beam. The three 
intensity profile monitors (IPMs), measured the intensity 
distribution of beam current in x and y. Two of the IPMs 
were coupled to a pair of fast ferrite cored steering mag- 
nets (FCSMs) that locked the beam path on the desired 
axis through the equipment. 
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FIG. 4: Method of scanning the CH2 blades through the 
beam. The blades measure 5 mm x 1.6 mm in cross sec- 
tion and extend 85 mm beyond the holders. Each PPM has 
two horizontally scanning and two vertically scanning blades. 
This figure shows a blade at the middle of a vertical scan. 

2. Intensity Profile Monitors 

The IPM signals were based on secondary emission 
from thin nickel foil strips arranged in harps placed be- 
tween 8 /jm thick aluminum foils (A very similar, ear- 
lier version is described in [24| and [25J). For IPM1 and 
IPM2 the nickel strips were 3 fj,m thick, 1.5 mm wide, 
and spaced 2.0 mm center to center. In IPM3, the nickel 
strips were 10 (im thick, 2.5 mm wide, and spaced 3.0 
mm center to center. Each IPM contained a vertical and 
a horizontal harp with 31 strips per harp. IPM1 also 
contained an aluminum normalization foil that provided 
a beam current signal to the liquid hydrogen target con- 
troller. The 31 signals from the strips of each of the six 
planes were individually amplified and digitized to pro- 
vide the beam intensity profiles in each spin state. Hard- 
ware beam centroid evaluators delivered signals propor- 
tional to the beam centroids. Corresponding correction 
signals were used to drive feedback loops to the pair of 
horizontal and vertical fast FCSMs (Fig|2J|. This allowed 
the beam intensity profile centroids to be kept within 
±l/im with a dc offset less than 50 fim from the desired 
axis. Sensitivities to helicity-correlated position and size 
modulations were determined with the beam unpolar- 
ized and with enhanced modulations introduced using 
the fast, FCSMs synchronized to the spin sequence. 



3. Polarization Profile Monitors 

The PPMs are described in detail in |2jj. Each PPM 
has four high density polyethylene (CH 2 ) blades, 5 mm 
x 1.6 mm in cross section. The blades were mounted 
on wheels rotating at 5 revolutions per second, and were 
scanned through the beam as shown in Fig. 0] two blades 
scanning vertically through the beam and two blades 



scanning horizontally. The figure shows one blade in 
the middle of a vertical scan. Protons elastically scat- 
tered from hydrogen in the blades were detected by sets 
of forward and recoil scintillator telescopes, the geome- 
try of which was set to select only pp elastic events at 
17.5° laboratory angle, near the maximum in the parity 
allowed analyzing power, A y . During vertical scans the 
up-down asymmetry was measured, and during horizon- 
tal scans the left-right asymmetry was measured. Infor- 
mation from each blade was stored in 80 time bins of 80 
/is each, for a total time window of 6.4 ms per blade. 
As the radius of rotation is 215 mm, this corresponds 
to ^(2-71") (215) mm = 43 mm centered on the nomi- 
nal beam axis. From the up-down or left-right asymme- 
try in each bin, a distribution of transverse polarization 
across the beam profile was generated. The two PPMs 
were operated with a 180° angular mismatch, producing 
8 equally spaced profiles in 1/5 second, or 40 profiles per 
second. (This is also the spin flip rate - see Fig. |3J) 

The master clock for sequencing the whole experiment, 
including hclicity changes, was derived from 2500-line op- 
tical encoders mounted on the PPM drive shafts. Each 
shaft encoder produced two 2500 cycle per revolution 
square waves in phase quadrature, permitting the direc- 
tion of rotation to be determined from whether one sig- 
nal was high or low at the positive going transition of 
the other. By recording both the rising and falling edges 
of both quadrature signals, the PPM control electronics 
obtained an effective resolution of 4 x 2500 = 10,000 
lines. The rotation of the two PPMs was synchronized 
by electronic gearing using position information provided 
by the shaft encoders. During steady state rotation at 5 
revolutions per second, the servo system was able to con- 
trol the angle between the two PPMs to ±3 milliradian, 
corresponding to one-third of the 1.6 mm target blade 
thickness at the location of the proton beam. 



4- Liquid Hydrogen Target 

The LH2 target had a flask 0.10 m in diameter and 
0.40 m long. The target scattered ~4% of the beam. It 
was important that the windows on the hydrogen vol- 
ume be flat and parallel to prevent motion of the proton 
beam from creating intensity noise. The hydrogen was 
contained by 25 /xm thick stainless steel windows. Up- 
stream and downstream of these inside the vacuum vessel 
were two helium filled chambers, 21 mm thick on the up- 
stream side and 10 mm thick on the downstream side, 
ensuring that the inner windows were flat. The outer 
windows of the helium cells were 50 (im thick stainless 
steel. The cryostat vacuum vessel was sealed with 25 fj,m 
thick copper windows. 

Although a thick target is desirable from the stand- 
point of increasing the signal, it was also very important 
that the tails of the beam profile be well contained in the 
150 mm x 150 mm aperture of the detectors; beam blow- 
up due to multiple Coulomb scattering limited the target 
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flask length, as well as the thickness of various entrance 
and exit windows and the thickness of the upstream IPM 
foils. With the target full, the beam size increased from 
approximately 4 mm (a of projected profile) at the center 
of the the upstream detector to 15 mm at the center of 
the downstream detector, and 22 mm by its exit. 

Rapid circulation (5 L/s) of the LH2 reduced density 
gradients and prevented boiling. A feedback loop us- 
ing fast and slow heaters controlled the target tempera- 
ture. The fast heater responded to a beam current signal 
from IPM1 (Sec. IIII D 2f> . and kept the heat load essen- 
tially constant when the beam current changed. The slow 
heater made fine adjustments to hold the LH2 tempera- 
ture at 19.3 K to within ±0.2 K over a several week data 
taking period. 



5. Transverse Field Ionization Chambers 

Each TRIC (a smaller earlier version is described in 
|27jl consisted of a cylindrical enclosure filled with 750 
liters of ultra-high purity hydrogen gas at a pressure of 
about 150 Torr, and contained an upper cathode plate 
operated at -8 kV, and a grounded lower signal plate. 
The TRICs also contained field shaping electrodes plus 
guard rings to ensure a uniform sense region, 150 mm 
wide by 150 mm high by 600 mm long between the par- 
allel electrodes. High-Z, low-energy spallation products 
from the chamber windows can cause large fluctuations in 
the signal, so the entrance and exit windows were located 
approximately 900 mm from the center to prevent spal- 
lation products from entering the active region. Other 
considerations in the design of the TRICs were noise due 
to delta ray (5-ray) production and ion pair recombi- 
nation. Recombination is reduced by lower gas pressure 
and higher voltage, both of which reduce the space charge 
density in the active volume. In practice, compromises 
must be made because too high a voltage causes increased 
noise from corona discharge and low pressure reduces the 
desired signal. 5-rays are electrons produced by collisions 
of the protons with the detector gas. The <5-ray signal is 
noisy, and the major contribution to this signal is from 
(5-rays at large angles. The transverse dimension of the 
ion chamber was a trade-off between a small transverse 
size which would minimize <5-ray noise and a large trans- 
verse dimension which accepts all the beam. The 150 
mm transverse dimension was the smallest that gave an 
acceptably low sensitivity to beam size modulation based 
on simulations. The main signal from the proton beam 
increases with the length of the active region. This was 
limited by the available space in the beamline; 600 mm 
was the longest practical length. 

The noise in the ion chamber signals has two incoherent 
contributions - shot noise from the statistical nature of 
the proton beam, and noise contributed by the chamber 
itself. The noise figure, a, expresses the chamber noise as 
a fraction of the shot noise. When the difference between 
the upstream and downstream chambers is taken, most of 



the shot noise contribution disappears because, except for 
the 4% scattered by the target, each upstream proton also 
passes through the downstream chamber. The chamber 
noise, on the other hand, does not cancel, and the run- 
ning time is dominated by chamber noise. In this experi- 
ment the counting time was approximately 15 times that 
which would have been expected from counting statis- 
tics alone. The fact that the run time was dominated 
by chamber noise resulted in the seemingly paradoxical 
result that better precision was obtained by lowering the 
beam current, because this reduced the detector noise fig- 
ure. In a series of test runs at progressively reduced beam 
currents, the A z distribution became narrower until ap- 
proximately 100 nA. The 200 nA selected for running was 
a compromise between lower current for better precision 
on A z and higher current for better precision from the 
PPMs, which were limited by counting statistics. 



IV. DATA ACQUISITION 

The data were taken in | second (200 ms) cycles, 
each cycle consisting of eight ^ second (25 ms) spin 

states arranged in the pattern (H 1 h H — ) or 

( — h H 1 h). This pattern makes the result in- 
sensitive to linear or quadratic drifts. The cycles were 
further arranged in an eight cycle (64 spin state) "super- 
cycle" , with the starting state of each cycle following the 

same (H 1 h +— ) or ( — h H 1 h) pattern. 

The initial spin state of each supercycle was chosen at 
random. In the frequency domain, this switching pattern 
contains odd multiples of 5 Hz, with the largest harmonic 
content at 15 Hz. The data acquisition produced on-line 
values for the amplitude and phase of the dominant 15 
Hz component. This was derived from a 16 spin state cy- 
cle. The 64 spin state supercycle gave the option of using 
more advanced digital filtering schemes on the difference 
signal, but these were only used during early develop- 
ment runs. 20% of the data were taken with all the spin 
flipping equipment running, but with the pumping lasers 
blocked with a shutter to guarantee zero polarization. In 
addition, as mentioned earlier, half of these "spin off" 
data were taken with an artificially enhanced intensity 
modulation synchronized with the spin flip. 

As shown in Fig. [3] each 25 ms spin state was di- 
vided into polarization measuring and asymmetry mea- 
suring intervals. Following a short dead time to allow 
the mechanical etalons which change the ion source spin 
state time to stabilize, was a 6.4 ms window during which 
one of the 8 blades of the scanning polarimcters passed 
through the beam. The TRIC and IPM signals were then 
integrated over exactly 1/60 second to eliminate sensitiv- 
ity to 60 Hz or harmonics of 60 Hz. The dead interval at 
the end of each state was intended as a buffer to absorb 
any timing jitter due to imperfect rotation speed of the 
PPMs. As it turned out, the timing jitter was less than 
0.1 ms and this buffer zone was more than adequate. 

The minimum data set for which a full set of helicity 
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25 ms (+ or - helicity) 



- 6.4 ms— 
PPM 



1/60 s 
TRIC 4c IPM 



FIG. 5: The timing of detector readout intervals within one 
spin state. The initial 0.8 ms allows for ion source settling 
time. The 6.4 ms PPM interval corresponds to a 43 mm 
scan centered on the beam axis. The TRIC interval was 1/60 
second to eliminate sensitivity to 60 Hz or harmonics of 60 
Hz. 



correlated beam properties could be extracted was a 0.4 
s "event pair" corresponding to two full 360° rotations 
of the PPMs, as this gave both spin states for each PPM 
blade. 

As noted in the ion source section, the spin state was 
transmitted as a frequency to prevent coupling the spin 
flip signal into the data acquisition. The effectiveness of 
this isolation was checked by running the complete data 
acquisition system including the ion source, but with 
detector signals supplied by a battery. False A z from 
electronic crosstalk was found to be less than 4 x 10~ 9 , 
an upper limit determined by the statistics of the test 
(crosstalk was probably less). 

The rotation speed of the PPMs is locked to a signal 
derived from the zero crossing of the 60 Hz AC line. Since 
the 0.2 s taken for one PPM rotation is exactly 12 cycles 
of the 60 Hz line, one would expect a given PPM blade to 
always pass through the beam at the same phase of the 
AC line. To prevent this, a small controlled phase slip 
was introduced. The rate of slip was programmable. For 
data taking it was set for one complete cycle of the 60 
Hz line in 18 minutes. Although the phase of one PPM 
relative to the other is "quantized" by the finite resolu- 
tion of the encoders, the angular frequency of rotation is 
not, and in principle can be locked to a square wave of 
any frequency. In practice, roundoff error in the motor 
control computer limited the choices of the reference fre- 
quency. The 18 minutes for 360° of slip corresponded to 
one of the acceptable values. 

In addition to the regular data taking runs, many ded- 
icated control measurements were made to measure the 
sensitivity of the apparatus to helicity correlated beam 
properties, specifically: position, size, intensity, trans- 
verse polarization, and energy modulation at the ion 
source. 



V. THE DATA 

A. Raw Data Set 

The parity data used in the final analysis were acquired 
during three major data taking periods - February- 
March 1997, July-August 1998, and May-June 1999. The 
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FIG. 6: E497 results before and after correction. The data 
are divided into 23 sets of alternating beamline tune. The top 
panel shows the results after beam quality cuts, but before 
the data were corrected for systematic errors. The bottom 
panel shows the data after correction. The error bars in the 
bottom panel are slightly larger due to uncertainties in the 
corrections. 

data were recorded as a series of "runs" of approximately 
one hour each. Statistics for the combined set are 

• 3.8 million event pairs 

• 375 runs in positive beamline helicity 

• 368 runs in negative beamline helicity 

• 80 position modulation runs 

• 81 size modulation runs 

• 40 energy modulation runs 

• 109 neutral axis scans 

The runs were grouped into 23 sets, a new set being 
started following each re-tune of the beamline (usually 
involving a reversal of the helicity at the parity apparatus 
relative to the ion source). 
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B. Reduced Data Set 

During data taking every effort was made to maintain 
the optimum beam conditions (quiet beam, good longitu- 
dinal polarization, low transverse polarization, low first 
moments of transverse polarization, beam on axis, stable 
position and size, low intensity modulation), and to stop 
data taking to correct the beam and cyclotron tune when 
these conditions were not met. Nevertheless, some data 
were recorded for which conditions were not ideal. The 
raw data sample was therefore subjected to data quality 
cuts prior to inclusion in the final analysis. 

The data quality cuts were conservatively chosen to 
ensure that corrections to the reduced data set based 
on critical hclicity-corrclatcd beam parameters would be 
consistent with sensitivities measured in ancillary cali- 
bration experiments. The data quality cuts first elimi- 
nated 46 data runs for which the TRIC difference signal 
was anomalously noisy, indicating unstable beam condi- 
tions. They also eliminated any data for which diagnostic 
monitor outputs could be considered spurious. The cuts 
used are summarized in Table ITU There were both "hard 
cuts" at a fixed value of beam parameter and "soft cuts" 
at ±3<7 from the mean value. The entire data reduction 
process reduced the size of the total set by 30%, but sig- 
nificantly improved the quality of the data sample. The 
xf, for the uncorrected data of the 23 sets went from 
11.3 per degree of freedom before the cuts were made 
to 5.5 after the cuts. It is important to emphasize that 
the approach used to determine the sensitivity to intrin- 
sic polarization moments, which was the dominant sys- 
tematic error correction, required consistent and stable 
beam conditions - in particular, the ratio of beam sizes 
at IPM1/IPM2 and the ratio of intrinsic polarization mo- 
ments at PPMI/PPM2 must be constant over the data 
sample used to determine the sensitivities, or the method 
is invalid. 

It is important to emphasize that no data cuts were 
applied to the parity violating asymmetry A z itself. The 
ultimate test of the consistency and validity of the data 
analysis and corrections procedure is the quality of the 
corrected data set shown in figure The top panel in 
Fig. shows the results for the reduced data set, before 
corrections were made for systematic errors. The bottom 
panel shows the corrected results from which the final 
A z was obtained. This corrected data set was subjected 
to a regression analysis in which residual sensitivities to 
hclicity correlated beam parameters were explored. No 
statistically significant residual sensitivities were found, 
lending confidence to the interpretation of the results. 
Details of the corrections procedure are given in the fol- 
lowing sections. 



VI. SOURCES OF ERROR 

Sources of error can be divided into those that are re- 
lated to the beam helicity and those that are not. Beam 



TABLE II: Summary of data reduction cuts. Except for the 
beam intensity, the center of the acceptance window for each 
parameter is the measured centroid of that beam parameter's 
distribution. The hard cuts are wide enough to include at 
least four standard deviations. 



Beam Parameter Acceptance Window 



Neutral axis x, y 

Position 

Width 

Skew 

Intensity 

Intrinsic Moments 



±0.3 mm 

±3cr (soft) 

±1.0 mm 
±0.2 mm 
196-204 nA 
±3.0 mm 



property changes that are synchronized with spin flip 
("coherent" or "helicity correlated" changes) can shift 
the centroid of the A z distribution and cause a false sig- 
nal of parity violation. Changes that arc not hclicity cor- 
related, such as detector noise and random variation in 
beam properties, do not bias the result, but they increase 
the run time required to reach a given precision. 



A. Random Changes 

The ultimate limit to the statistical precision of the 
experiment was set by the counting statistics of the scat- 
tered protons. For a target scattering 4% of the 200 nA 
beam, the rate of scattered protons is 50 GHz, which 
would make it possible to measure A z to ±0.2 x 10 -7 in 
20 hours if individual scattered protons could be counted. 
However, 50 GHz was too high for direct counting, and so 
current mode detection was used. As noted already in the 
ion chamber section, the counting time was dominated by 
detector noise. Other random variations in beam proper- 
ties such as intensity and position also contributed noise 
and further increased the required run time. It proved 
to be a net advantage to devote significant beam devel- 
opment time to producing quiet, stable, beam. 



B. Helicity Correlated Changes 

The approach to minimizing and correcting for the 
false A z signal due to helicity correlated beam property 
changes was as follows: 

• Careful design and operation of the TRIUMF op- 
tically pumped polarized ion source and cyclotron 
made it possible to change the spin direction with 
very little effect on the other beam properties. 

• The design of the parity equipment and the oper- 
ating conditions of the experiment were carefully 
chosen to minimize the sensitivity to helicity corre- 
lated changes. 
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• Calibration runs determined the sensitivity to hc- 
licity correlated modulations. 

• The beam properties were continuously monitored 
during data taking so the actual helicity correlated 
changes were known, and appropriate corrections 
were applied. 



VII. CORRECTION FOR SYSTEMATIC 
ERRORS 

The effects of modulations in beam intensity, position, 
size, transverse polarization, and energy were considered. 

Measurements made early in the experimental program 
revealed significant correlations between beam parame- 
ters. For example, horizontal and vertical beam motion 
are often highly correlated. To be sure of extracting the 
correct sensitivities, separate control measurements were 
made in which each beam property in turn was artificially 
modulated, and the effect on apparent A z was recorded. 
The sensitivities obtained by this method were consistent 
with sensitivities extracted by multilinear regression us- 
ing the natural variation of beam parameters. Most im- 
portantly, the regression analysis showed no significant 
sensitivity to products of beam modulations, so the false 
A z was the sum of the contributions of the individual 
coherent modulations. 

For the sensitivity calibrations made by varying one 
parameter at a time, measurements were made using 
modulations of different size and sign. The false asym- 
metry was linear in the modulation, passing through zero 
at zero modulation. To determine the sensitivities used 
for data reduction, large values of modulation were used 
to give accurate results in a short time. To determine 
the sensitivity to first moments of transverse polariza- 
tion, which gave the largest false asymmetry, much larger 
transverse polarizations were used for calibration than 
were present in the parity beam. This was justifiable be- 
cause false asymmetries must be linear in polarization, 
and the use of large polarizations allowed the displace- 
ments needed for the calibration to be reduced. The mea- 
sured false asymmetry was linear in polarization moment. 



A. Coherent intensity modulation (CIM) 

As the energy of the beam was 27 MeV lower at TRIC2 
than at TRIC1, the TRIC2 gas pressure was lowered 
relative to TRIC1 to equalize ion chamber currents. A 
hardware gain on the TRIC1 signal was then adjusted to 
minimize common-mode noise in the amplified difference 
signal, 1000 x (Ja — gl\) which was digitized, where g is 
the hardware gain; g is very close to unity. The helicity 
correlated part of the difference, which should be pro- 
portional to A z , nevertheless contains a coherent current 
modulation component due to imperfect common mode 
rejection arising from gain mismatch and nonlinearity of 



the ion chambers. If the current signal from TRIC2 is 
I2 = al^(l — 5(1 ± P Z A Z )) then, expressing the differ- 
ence in current for the two helicity states as 7 X = I\ ±5h , 
one can define an "analog asymmetry" , 

_ (J+ - git) - (I 2 - g/f) 

= SaP z A z + (aT g + hT^j ^, (1) 

where S is the nuclear scattering probability in the LH 2 
target (S ~ 0.04) and T = 1 — S. The quantity a is a 
function of TRIC1 and TRIC2 gas gains, and is nominally 
equal to 1/T; it was adjusted with TRIC2 gas pressure. 
The hardware gain g was set by zeroing the false parity 
signal when the beam current was modulated by the pho- 
todetachment laser. Because of the nonlinearity term, g 
must be reset if the beam current changes. 

A beam current stabilization system was implemented 
which restricted beam current excursions to the range 
200 ± 2 nA 01 ■ The current stabilization system took 
a current signal from the upstream ionization chamber 
(TRIC1) and fed a correction voltage to an electrostatic 
quadrupole lens just upstream of a set of slits in the in- 
jection line. The injection line was tuned with a slight 
excess of current so that, at the set point, about 10% of 
the beam was skimmed by the slits. To exclude coupling 
of current modulation to position and energy modula- 
tions, the sampling rate was slow (0.5 Hz), much slower 
than the spin- flip rate of 40 Hz. The active range of 
the current stabilization loop was quite low (±10%) to 
avoid correcting drifts caused by large excursions of cy- 
clotron tune. The current loop operation was a sensitive 
indicator of cyclotron stability, and operators made ad- 
justments to the cyclotron tune when necessary. Even 
under conditions of constant beam current, periodic ad- 
justments of g were needed to compensate for small drifts 
in the ion chamber gains. The setting for best common 
mode rejection was checked periodically by turning on 
the photodetachment laser, creating a (~ 0.1%) coherent 
intensity modulation, and performing a scan of subtrac- 
ter box settings. The results of such a calibration are 
shown in Fig. [7] Because the difference signal was multi- 
plied by 1000, the DC level was a sensitive indication of 
the setting of the hardware fine gain. A fractional change 
of 10~ 4 in the hardware fine gain caused a DC shift of 
0.7 V in the subtracter box output. 

During data taking, coherent intensity modulation 
(CIM) was measured continuously by the ion chamber 
(TRIC1) upstream of the target. As noted in section 
IIII Bl this was normally less than a few parts in 10 5 . Be- 
cause the sensitivity to CIM changed between the dedi- 
cated CIM runs used to re-set the subtractor box, the sen- 
sitivity to CIM, like the CIM itself, was measured during 
data taking. Periods of enhanced ( 0.1%) coherent inten- 
sity modulation were interleaved with the main data. A 
typical sensitivity was A{ alse = (1.6 x 10~ 4 )(/ b -200)(^) 
where is the beam current in nA. Table ITTT1 shows the 
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FIG. 7: A scan of the subtractor box settings taken during 
the July 1998 running period. The plot shows the sensitivity 
to CIM as a function of the subtractor box DC level. The DC 
level is a sensitive indication of the setting of the fine gain 
control; a fractional gain change of 10 -4 causes a shift of 0.7 
V in the DC level. This particular scan gave an ideal DC level 
of -1.9 V for zero sensitivity to CIM. 
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FIG. 8: CIM data for a typical one hour run. The plot shows 
a histogram of the sensitivity to . Data have been excluded 
for which |^| < 0.06%. The mean is 0.23 and the standard 
deviation 2.6, in units of (10 _4 /%). 



average sensitivities to coherent intensity modulation for 
the different running periods, and Fig. |H| shows a his- 
togram of the sensitivities for a typical (~ 1 hour) run. 



B. Coherent position modulation 



TABLE III: Average sensitivities to Coherent Intensity Mod- 
ulation. The sensitivities are determined from all CIM data 
with > 0.06 %. 



Set 


Sensitivity (10 


7%) CIM events 


1997 


36 ± 19 


83,636 


1998 


138 ±8 


106,520 


1999 


73 ±11 


90,955 



and 1999 data, IPM1 and IPM2, located 1.8 m apart 
along the beamline, were used. Using two IPMs permits 
measuring both tilts and parallel shifts of the beam. The 
measured coherent modulations are shown in Tables IIVI 
and0 The sensitivity to beam motion was measured in 
separate control measurements during which the beam 
position was modulated in a variety of ways using the 
FCSMs. 

The false analyzing power arising from helicity corre- 
lated beam position was parameterized as 



AA Z = 



Ax i + Ax2 



Axi — Ax2 



A Vl + Ay 2 



^ Ayi - Ay 2 



(a x xi + b x x 2 + Cx) 
(d x xi + e x x 2 + f x ) 
(ayyi + b y y 2 + c y ) 
(dyij! + e y y 2 + f y ) 



(2) 



where Ax and Ay are the horizontal and vertical helicity 
correlated beam motion, x and y are the beam position 
and the 1 and 2 subscripts refer to IPM1 and IPM2. The 
parameters a to f were extracted from a fit to the cali- 
bration data. Details may be found in Ref. Figure 
shows the results of such a fit. The abscissa shows the 
false asymmetry predicted using Eq. 2 with the fitted 
parameters, and the ordinate shows the false asymmetry 
actually measured with the calibration data. 



C. Coherent beam size modulation 

Coherent beam size modulation was also measured by 
the two IPMs, and the observed values are listed in Ta- 
bles ^3 and[V] The sensitivity to beam size was deter- 
mined by control measurements in which the beam size 
was intentionally modulated by driving the FCSMs as 
quadrupoles. 

The false analyzing power due to beam size modulation 
was expressed as 



Coherent position modulation was measured by two of 
the three IPMs. During the 1997 running period, IPM1 
and IPM3 were used because their separation was greater 
and, all else being equal, position control should have 
been better. Unfortunately IPM3 picked up noise from 
the liquid hydrogen target circulation fan, so for the 1998 



AA Z = a x a Xl Aa Xl + /3 x a X2 Aa X2 

+ a y <Ty 1 A(Jy 1 + f3y(Ty 2 A(Jy 2 (3) 

where a is the RMS beam size and Ac is the helicity 
correlated change in beam size. As with position mod- 
ulation, the parameters were extracted from a fit to the 



11 



1.0 



0.5 - 



< 



0.0 



-0.5 - 



-1.0 



I 




J*' 





-1.0 -0.5 0.0 0.5 1.0 

Fitted Value (10 ) 

FIG. 9: Calibration of sensitivity to coherent position mod- 
ulation. The abscissa shows the false asymmetry based on 
sensitivities extracted from calibration runs and the ordinate 
shows the asymmetry actually measured. 




TABLE IV: Summary of Position and Size Modulation for 
the 1997 data set. These data used IPM3 for fast position 
control. Changing the Wien filter sets the spin direction in 
the cyclotron relative to the ion source and the beamline sets 
the helicity at the parity apparatus relative to the cyclotron. 





Wien 


filter ± 


Wien filter — 


Beam Line: 


+ 




+ 




Aa xi (nm) 


119 ±49 


-137 ±43 


-101 ±59 


8 ±48 


Aa X3 (nm) 


65 ±79 


-156 ±76 


-110 ± 118 


-45 ±91 


Aa yi (nm) 


126 ±44 


-212 ± 44 


-36 ± 54 


41 ±45 


Aa V3 (nm) 


81 ±75 


-200 ± 72 


-54 ± 116 


-11 ±84 


Axi(nm) 


5 ±33 


19 ±28 


-34 ± 40 


38 ±36 


Ax$(nm) 


-30 ± 34 


-40 ± 34 


10 ±46 


4 ±37 


Ayi(nm) 


-46 ± 26 


89 ±27 


-7 ±33 


-10 ±27 


Ay s (nm) 


-27 ± 36 


-12 ±33 


24 ±55 


-78 ± 43 



TABLE V: Summary of Position and Size modulation for the 
1998 and 1999 data sets. IPM2 was used for fast position 
feedback, replacing IPM3. The Wien filter was ± for all these 
data, producing spin-up in the cyclotron. The beamline ro- 
tates this to ± or — helicity at the parity apparatus. 



Beam Line: 



± 



Aa xi (nm) 
Aa X2 (nm) 
Aa yi (nm) 
A(j V2 (nm) 
Ax\(nm) 
Ax2(nm) 
Ay! (nm) 
Ay2(nm) 



-17 ±38 

8±8 
-6 ±37 
17 ±7 
1 ±14 
7 ± 11 
17 ±11 
16 ±10 



0±35 
-35 ±7 
1±35 
-37 ±5 
-1±6 
-9 ±5 
3±8 
2± 7 



'-0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 
Fitted Value (10 4 ) 

FIG. 10: Calibration of sensitivity to coherent size modu- 
lation. The abscissa shows the false asymmetry based on 
sensitivities extracted from calibration runs and the ordinate 
shows the asymmetry actually measured. 

calibration data. Because the quadrupolc configuration 
of the FCSMs caused the beam to be steered somewhat 
when the size was modulated, correction for position 
modulation had to be included in the fit. Figurc lTUl shows 
the results of a fit to a series of size modulation calibra- 
tion runs. 



D. Transverse polarization 

If the proton spin is not perfectly longitudinal, the 
small transverse component will reverse with helicity. 
This can couple to the relatively large (~ 0.3) parity 
allowed analyzing power to cause a false parity violating 
signal. Figurc lTTl uscs a pencil beam with a vertical trans- 
verse polarization to illustrate the mechanism. Because 
the parity allowed analyzing power is positive, slightly 



more beam will be scattered to the left in the positive 
helicity state, and to the right in the negative helicity 
state. If the beam passes through the center (neutral 
axis) of the detector, as shown in the left-hand panel, 
the response will be the same for both helicity states and 
no false effect arises. If, on the other hand, as shown in 
the right-hand panel, the beam does not pass through 
the center of the detector, more signal will be recorded 
in one helicity state than in the other. To a very good 
approximation, the effect is found to be proportional to 
the size of the transverse component multiplied by the 
distance the beam is off center at the detector - i.e. to 
the first moment (x)(P y ) of transverse polarization at the 
detector. 

In a field-free region, a real particle beam of finite ex- 
tent is made up of a bundle of straight rays like the pencil 
beam in this example. The first moment for the beam is 
the linear sum of the first moments of the individual rays. 
Since the first moment for a ray is proportional to the 
distance of the ray from the zero axis and varies linearly 
with distance along the beam, changing sign where the 
ray crosses zero, the first moment for a particle beam will 
vary linearly with distance along the beam. A real beam 
can have a substantial first moment of transverse polar- 
ization even if the net transverse polarization is zero. In 
the example of Fig. ^2 the vertical polarization could be 



12 




detector sensitive 
volume 



scatterer 




detector sensitive 
volume 



scatterer 



FIG. 11: The effect of transverse components of polarization. 
The combination of a transverse polarization component with 
a beam displacement from the neutral axis will cause a false 
signal of parity violation. 
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FIG. 12: A horizontal neutral axis scan. A vertically polarized 
beam (P y = 0.76) was scanned horizontally and the sensitivity 
to P y plotted as a function of x projected to the "magic z" 
[z location where the false A z does not depend on angle). A 
beam with x pro j — —2.7 mm was found to produce no false 
signal from vertical polarization components. 



"up" on the right side of the beam and "down" on the 
left side. 

To determine the first moment sensitivities, test runs 
were made with pure vertical and pure horizontal po- 
larization. By scanning the vertically polarized beam 
horizontally it was possible to generate known moments, 
(x)(P y ) 1 and, by scanning the horizontally polarized 
beam vertically, known (y)(P x ) moments could be gen- 
erated. Although in principle the first moments must be 
known at the two detectors, because first moments vary 
linearly with distance along the beamline, knowledge of 
the first moments at any two locations is sufficient. 

The transverse polarization profiles were measured by 
the two Polarization Profile Monitors (PPMs in Fig. |5J) 
located 1.8 m apart along the beamline and, for conve- 



nience, the first moment sensitivities were expressed in 
terms of the sensitivities at the PPMs. The first mo- 
ment scans also defined the beam trajectory for which 
there was zero sensitivity to average transverse polariza- 
tion. Such a trajectory was referred to as a polarization 
neutral axis. Polarization neutral axes were determined 
in both the horizontal and vertical directions and, dur- 
ing data taking, the beam was held on this neutral axis 
by the servo system. Figure ^] shows an example of a 
horizontal neutral axis scan. Another consequence of the 
linear behavior of the first moments with distance, z, is 
the presence of a "magic z" downstream of TRIC2. Any 
beam with zero first moment at this z will cause no false 
effects. Also notice that the measurements shown in Fig. 
1121 confirm the linear relation between first moment and 
false A z . 

In analyzing the effect of first moments, extrinsic first 
moments, (x)(P v ) and (y)(P x ), caused by a beam which 
has some net transverse polarization and whose centroid 
is displaced from the neutral axis, were treated separately 
from intrinsic first moments, (xP y ) and (yP x ), which do 
not depend on the position of the beam, but rather arise 
from the distribution of components within the beam. 
By holding the beam on the polarization neutral axis it 
was possible to virtually eliminate corrections for extrin- 
sic first moments. Intrinsic first moments, on the other 
hand, are independent of beam position, arise in the cy- 
clotron and beamline, and were very hard to control. 

The false asymmetry arising from extrinsic and intrin- 
sic first moments can be written 



AA, 



(P x )[ ai {( yi )-yr) +a 2 {(v2) 



(p v )M(xi) - 

+ai(yP x )i 
+b 1 (xP y ) 1 



...i la 
1 



) + 



b 2 ((x 2 ) - 
+a 2 {yPx)2 



-botxR 



y)2 



(4) 



where the first two lines are the contribution from ex- 
trinsic first moments, and the second two lines are the 
contribution from intrinsic first moments. Note that the 
first moment sensitivities, ai, a 2 , b\, and b 2 , are identical 
for intrinsic and extrinsic moments. 

The results of the neutral axis scans are shown in Table 
EH The quantities (-a iy ™ - a 2 y^ a ) and (-bix^ a - 
62^2°) depend on what axis x and y are measured from. 
They are zero if x\, yi, x 2 , and y 2 are measured relative 
to the neutral axis. One notes that the sensitivities to 
first moments are very consistent from run period to run 
period. In units of 10 -7 per ^m, the sensitivity to (y) (P x ) 
is 1.8 at PPM1 and -2.5 at PPM2. The sensitivity to 
(x)(P y ) is -1.5 at PPM1 and 2.0 at PPM2. 

There were no magnetic elements after the first PPM, 
so the first moments of transverse polarization varied lin- 
early with position along the beamline. Furthermore, for 
a fixed setting of the upstream beamline magnets, the 
first moments at PPM1 and PPM2 tended to scale to- 
gether so that, over a wide range of first moments, the 
ratio of first moments at PPM1 and PPM2 had a con- 
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TABLE VI: Fitted first moment sensitivities from the neutral 
axis scans. 



Sensitivity (1CT 7 ) 


Feb97 


Jul98 


May99 


oi(fim _i ) 


1.91 ±0.03 


1.85 ±0.01 


1.83 ±0.02 




-2.47 ±0.01 


-2.47 ±0.01 


-2.44 ±0.01 


(-a lV r - a 2V r) 


-0.55 ±0.01 


-1.15 ± 0.01 


-0.26 ± 0.01 


h(ij,m 1 ) 


-1.51 ±0.03 


-1.49 ± 0.01 


-1.49 ± 0.02 


&2(M m_1 ) 


2.02 ±0.04 


2.00 ±0.01 


1.98 ±0.01 


(— bix 1 — b 2 x 2 ) 


1.31 ±0.03 


-0.35 ± 0.01 


0.50 ± 0.01 


xl 


1.17 


1.96 


2.05 



TABLE VII: Settings of the Wien filter and beamline for the 
23 data sets. 



set, 


Wien filter, beamline) 




1,(±, 


+ 


9,(+,-) 


17, (+, 


i ) 


2,(±, 




io,(±,±) 


18, (+, 




3,(±, 


+ 




19, (+, 


±) 


4,(+, 




12,(±,±) 


20, (+, 




5,(+, 


+ 


13, (+,-) 


21, (+, 


±) 


6,(- 




H (+,+) 


22, (+, 




7,(- 




15, (+,-) 


23, (±, 


+) 


8,(±, 


+ 


16,(±,±) 







stant value. Since the neutral axis scans showed that the 
first moment sensitivity at PPM1 was of opposite sign 
to, and 75% of the magnitude of, that at PPM2, if the 
PPM2 first moment was always 75% of the PPM1 first 
moment, then the effective sensitivity to first moments 
would be zero. This would correspond to a beam whose 
first moment of transverse polarization goes through zero 
at the magic z. By adjusting the beam convergence, an 
attempt was made to achieve this magic moment ratio 
between PPM1 and PPM2. Unfortunately, the ability of 
the PPMs to measure this ratio to sufficient precision in 
a reasonable amount of time was very limited, so the can- 
cellation was not perfect and the residual, effective, first 
moment sensitivity had to be extracted from the data. 
This could be done cither by regressing the first moment 
ratio from the data and using the sensitivities measured 
in the neutral axis scans, or by regressing the correlation 
of A z with average first moment directly from the data. 
The two methods agreed with each other, but the latter 
method produced the smallest statistical spread and was 
the method used to correct the data. 



E. Energy Modulation 

Since the beam energy in the downstream TRIC was 
on average 27 MeV lower than in the upstream TRIC 
due to energy loss in the target, hclicity-correlated en- 
ergy modulation caused a systematic error due to the 
nonlinear energy dependence of the proton beam energy 
loss in the hydrogen gas of the TRICs. The sensitiv- 
ity to coherent energy modulations was determined us- 
ing an RF accelerating cavity placed upstream of IPM1 
in the beam line. The measured sensitivity of false A Zl 
(2.9 ±0.3) x 10~ s eV~ 1 , was in excellent agreement with 
predictions based on the variation of stopping power with 
energy. 

Energy modulation of the extracted beam was caused 
by position modulations of the radial intensity distribu- 
tion at the stripping foil; this converted radial position 
modulation of the injected beam to energy modulation of 
the extracted beam. The primary coherent energy mod- 
ulation produced in the source was converted to position 
modulation in the injection beamline and then back to 
energy modulation at the extraction foil. Direct mea- 



surements using a magnetic spectrometer (1.2 m disper- 
sion) in another beamline (4B) at TRIUMF showed the 
energy modulation of the extracted beam to be approxi- 
mately 100 to 200 times greater than the energy mod- 
ulation at OPPIS. During the parity runs this direct 
measure could not be made, but frequent measures of 
energy modulation at OPPIS and of dA z /dE opp i S were 
made. The dA z /dE opp i S sensitivity was measured by 
applying a square wave voltage of 0.5 V amplitude to 
the electrically isolated sodium ionizer in OPPIS. These 
dA z /dE opp i S runs were then corrected for all known sys- 
tematic errors, (dominated by dl/I), and it was assumed 
that the residual false A z arose from energy modula- 
tion of the extracted beam. A comparison of the mea- 
sured dA z /dE opp i S to the dA z /dE measured in the beam- 
line with the RF cavity (Sec. IV11 Ffl indicated that the 
cyclotron amplified dE opp i S by a factor of about 130, 
in agreement with the magnetic spectrometer measure- 
ments. 

The primary energy modulation caused by the optical 
pumping lasers was measured by using an electrostatic 
beam energy analyzer in the polarized source and an in- 
tensity profile monitor with 16 collector strips 2.5 mm 
wide, and with 3.0 mm spacing to measure beam po- 
sition modulation downstream of the steering analyzing 
plates. The monitor was mounted on a remotely con- 
trolled swinging arm. Two measurements of coherent 
position modulation for the right and left monitor posi- 
tions allowed separation of the energy and position mod- 
ulation components of the OPPIS beam. An accuracy of 
0.2 meV could be achieved in ten minutes of integration 
time. The modulation magnitudes were quite sensitive 
to the pumping laser asymmetry between the two po- 
larization states; after careful laser tuning, the coherent 
energy modulation was reduced to 1-2 meV and the co- 
herent position modulation to the 20 nm level. 

Although the frequent measurement of AE and 
dA z j 'dE ppi S helped to set limits on the false A z from 
energy modulation of the extracted proton beam, it is 
significant that the energy modulation of the extracted 
proton beam could not be measured directly at the parity 
apparatus. To cancel its effects, use was made of the fact 
that when the beamline hclicity is reversed, the effects of 
true A z reverse, but the effects of energy modulation do 
not. This is because energy modulation can only arise 
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TABLE VIII: Overall corrections for systematic errors. The 
table shows the average value of each coherent modulation, 
the net correction made for this modulation, and the uncer- 
tainty resulting from applying the correction. 



Property 


Average Value 


10 7 AA Z 






(correction) 




1.68 ± 0.29(stat.) 




y * P x (fj,m) 


-0.1 ±0.0 


-0.01 ±0.01 


x * P,,(am) 
y 1 


-0.1 ±0.0 


0.01 ±0.03 


(yP x )(um) 


1.1 ±0.4 


0.11 ±0.01 


{xP y ){um) 


-2.1 ± 0.4 


0.54 ± 0.06 


AI/I(ppm) 


15 ± 1 


0.19 ±0.02 


position + size 




0±0.10 


AE{meVat OPPIS) 


7-15 


0.0 ±0.12 


electronic crosstalk 




0.0 ± 0.04 


Total 




0.84±0.17(syst.) 


A c ° rr (lQ-') 
X l(23sets) 


0.84 ± 0.29{stat.) ± 0.17(syst.) 
1.08 



in the ion source and cyclotron. The magnets used to 
rotate the spin from "up" to positive helicity or "up" to 
negative helicity arc downstream of all sources of energy 
modulation and do not affect the beam energy. Data 
were generally taken with alternating beamline helicity, 
and during the 1997 run data were also taken with the 
Wien filter reversed, an independent method of reversing 
the proton beam helicity relative to the ion source. The 
Wien filter and beamline settings are detailed in Table 
IVIII That the xt f° r the 23 sets is only 1.08 following cor- 
rections, shows that the effects of uncorrected systematic 
errors, including energy modulation, were small. Care 
was taken to balance the amount of data taken in the 
two beamline helicities so that if some false A z from en- 
ergy modulation was present, it would tend to cancel in 
the final weighted average. Information on the energy 
modulation at the ion source and the sensitivity to this 
modulation was used to include the effects of uncorrected 
energy modulation in the error budget, as described in 
more detail in the next section. 



VIII. METHOD OF APPLYING THE 
CORRECTIONS 

The false A z arising from a given coherent modulation 
was found by multiplying the measured modulation by 
the sensitivity to that modulation. Uncertainties from 
the variance of the corrected A z distribution are referred 
to as "statistical" . The uncertainty quoted is the stan- 
dard error in the mean of the A z distribution. Uncertain- 
ties in A z resulting from uncertainties in the sensitivities 
are referred to as "systematic" because an incorrect sen- 
sitivity will cause a systematic shift in the mean of the 
A z distribution. These errors are, however, statistical 
in nature, as they arose from statistical uncertainties in 



the knowledge of the sensitivities to the various coherent 
modulations. 

Table I Villi summarizes the overall corrections to the 
parity data. To produce the 23 corrected A z distributions 
shown in Fig. H3 the following procedure was followed. 

1. The data in each set were grouped into bundles of 
10000 event pairs per bundle. 

2. Each bundle was corrected according to the ob- 
served coherent modulations (except position and 
size) for that bundle, giving a corrected A z for each 
bundle. The variance of corrected A z values in a 
set determined an error bar for that set. This is re- 
ported as the "statistical" uncertainty. No correc- 
tions were made for position and size, because the 
net correction was consistent with zero and when 
the corrections were applied it was found that they 
slightly increased the residual correlations of A z 
with position and size, as well as increasing the 
variance of the corrected A z distribution. 

3. The uncertainties in the nil correction for position 
and size modulation were added in quadrature to 
each of the 23 sets. This is included in what is 
reported in Table fVTTTl as "systematic" uncertainty. 

4. The uncertainties resulting from uncertainties in 
the various sensitivities were added in quadrature 
to the corresponding data. The uncertainties in the 
sensitivities are independent of what is accounted 
for in step 2, so it is justified to add them in quadra- 
ture to obtain the total error bar on each of the cor- 
rected A z for the 23 sets. This uncertainty is also 
included in the "systematic" uncertainty in Table 

Enu 

5. The A z reported is the weighted mean of the 23 
data sets with a weight 1/err 2 , where err is the 
"total" uncertainty, not including energy modula- 
tion. The error bars shown in the bottom panel of 
Fig. EJare these "total" uncertainties. 

6. As mentioned earlier, the correction for extracted 
energy modulation was complicated by the fact 
that no direct measurement of energy modulation 
could be made in the parity beamline. The energy 
modulation sensitivities depended on the beamline 
tune, so the distribution of measured OPPIS en- 
ergy modulation values and dA z / dE opp i S sensitivi- 
ties was examined for a given beamline tunc, and 
a correction and a "worst case" uncertainty in the 
correction were estimated. 

Net corrections were calculated for energy modula- 
tion on a year by year basis. 1998 and 1999 required 
two corrections each, one for each beamline helicity. 
1997 needed four corrections, as two Wien Filter 
settings were used. Finally all the energy modula- 
tion corrections were combined and one correction 
was applied to the final A z . 
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FIG. 13: Theoretical predictions for A z and the most precise 
experimental data at 13.6 MeV (Bonn), 45 MeV (PSI) and 221 
MeV (TRIUMF). The solid curve shows the results obtained 
by Carlson et al. by adjusting the weak coupling constants 
for the best fit to the experimental data. 

The energy modulation correction shown in Table 
I Villi is the net effect of energy modulation over the 
three runs. The net correction was zero. The ±0.12 
uncertainty comes from the quadrature sum of the 
uncertainty in energy modulation plus the uncer- 
tainty in the energy modulation sensitivities. This 
was not included in the individual error bars for the 
23 A z numbers because the energy modulation was 
not known well enough. As a result, the reduced 
chi squared of 1.08 for the 23 sets is larger than it 
really "should" be if the energy modulation uncer- 
tainty was determined individually for each of the 
23 sets and included in the individual set by set 
error bars. 



IX. RESULTS 

The overall result for A z from the 23 sets is sum- 
marized in Table IVIIII After correcting for systematic 
errors, the longitudinal analyzing power is found to be 
A z = (0.84 ± 0.29(stai.) ± 0A7(syst.)) x 10~ 7 at 221.3 
MeV incident proton energy and the target and detector 
geometry of this experiment. Correcting for finite geome- 
try and target thickness (see Sec. IIII C 21 and combining 
the errors in quadrature, gives A z = 0.86 ± 0.35 at 225 
MeV for comparison with theoretical calculations. 

Parity violation in pp scattering has already attracted 
considerable theoretical interest, and many calculations 
of A z have been made [H El El El El These calcu- 
lations are shown in Fig. Altogether with the TRIUMF 
result (corrected to the 1 So — 3 Po zero crossing energy) 
and the most precise results at 13.6 MeV @ and 45 MeV 
0. Theoretical predictions for A z from several models 
are shown. The model of Driscoll and Miller [2^ is based 




FIG. 14: Present constraints on the weak meson- nucleon cou- 
plings based on the experimental data and recent calculations 
by Carlson et al.|28||. The bands are the constraints imposed 
by different experiments (Bonn 13.6 MeV, dashed; PSI 45 
MeV, solid; TRIUMF 221 MeV, dotted). The filled square 
and dotted rectangle are the DDH "best guess" and "reason- 
able range" respectively. Also shown are the 68% and 90% 
C.L. contours. 



on the Bonn potential to represent the strong N-N inter- 
action, together with the weak meson-nucleon coupling 
constants as given by Desplanques, Donoghuc, and Hol- 
stein (DDH) 0]. The prediction of Iqbal and Niskanen 
E3 has a A isobar contribution added to the Driscoll 
and Miller model on a semi-ad-hoc basis. The theoretical 
prediction of Driscoll and Meissner [3l| is based on a self- 
consistent calculation, with both weak and strong vertex 
functions obtained with a chiral soliton model. Finally, 
the quark model calculation of Grach and Shmatikov |32( 
takes explicit account of quark degrees of freedom. None 
of these predictions are in good agreement with the data, 
although they all have similar shapes due to the energy 
dependence of the strong interaction. 

A major source of uncertainty in these calculations is 
the value of the weak meson-nucleon couplings. Starting 
with the work of Desplanques, Donoghue, and Holstcin 
(DDH) ! 2| , many theoretical calculations of these quan- 
tities were made El |3 El Ell but theoretical uncer- 
tainties remained large. Until the TRIUMF experiment, 
the coupling constants were very poorly constrained by 
experiment. 

Figure ITU shows the limits on the weak meson-nucleon 
couplings hP p and h v J now imposed by the low energy 
results 0, H and the present TRIUMF result. The er- 
ror bands are based on a calculation by Carlson et al. 
[28| assuming the Argonne vir (AV-18) potential |37| . 
the Bonn 2000 (CD-Bonn) [38| strong interaction cou- 
pling constants, and including all partial waves up to 
J=8. Although the TRIUMF measurement is not sensi- 
tive to A z from SP mixing, and the contribution from 
PD mixing contains no h^f contribution, there is some 
h^f dependence arising from the higher partial wave mix- 
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ings. The net result is that the acceptable band de- 
fined by the TRIUMF measurement is almost orthogo- 
nal to that defined by the low energy measurements, and 
greatly reduces the acceptable ranges of both hf and 
h v J . Adjusting these coupling constants for the best fit 
to the pp data, including the TRIUMF 221 MeV point, 
Carlson et al. H| estimate hf = -22.3 x 1(T 7 and 
hPJ = 5.17 x 10" 7 , compared to the DDH "best guess" 
values of hf = -15.5 x 10" 7 and hPJ = -3.0 x 1(T 7 . The 
solid curve in Fig. El is calculated using the Carlson et 
al. adjusted couplings. 

The reduction in the experimentally allowed range for 
weak meson-nucleon coupling constants also has impli- 
cations for the analysis of electroweak radiative correc- 
tions in backward angle parity violating electron scatter- 
ing. By combining back angle electron scattering data for 
hydrogen ^ and deuterium 0, the SAMPLE collabo- 
ration was able to extract values for both the isovector 
axial e-N form factor, G e A {T =1), and the strange mag- 
netic form factor, G S M at a momentum transfer Q 2 = 0.1 
(GeV/c) 2 . Their value G e A (T = 1) = +0.22 ± 0.45 ± 0.39 
differs significantly from the value G e A (T = 1) = — 0.83± 
0.26 arrived at by Zhu et al. Q by applying one quark 
and many quark (proton anapolc moment) radiative cor- 
rections to the well known nucleon axial charge as mea- 
sured in neutron /3-decay. hf enters in the calculation 
of the proton anapole moment, but its contribution is 
small and one would require hf ~ —180 x 10~ 7 to bring 



the value up to the +0.22 of the SAMPLE measurement. 
Such a value of hf is now clearly ruled out by the pp 
parity violation data. 



X. CONCLUSION 

The parity violating analyzing power, A z , in pp elas- 
tic scattering has been measured at 221.3 MeV incident 
proton energy. The result constrains theoretical calcu- 
lations of A z in an energy region not previously covered 
experimentally. In the case of meson exchange calcula- 
tions, it constrains principally the value of the weak rho 
meson- nucleon coupling hf, but, when the low energy pp 
data are included, strong constraints are placed on the 
acceptable values of both hf and K^f . This result has 
implications for the interpretation of other experiments. 
For example, it rules out incorrect values of these cou- 
plings as an explanation for the disagreement between 
the SAMPLE isovector axial form factor result and 
the calculation of Zhu et al. . 
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